The Evolution of the Pion Distribution Amplitude 
in Next-to-Leading Order 



D. Midler 
Stanford Linear Accelerator Center, 

■^f- . Stanford University, Stanford, California 94309 

On 
On 

> 

O ■ Abstract 



ON 



> 

oo 



The evolution of the pion distribution amplitude in next-to-leading order is studied 
for a fixed and a running coupling constant. In both cases, the evolution provides a 
logarithmic modification in the endpoint region. Assuming a simple parameterization 
of the distribution amplitude at a scale of Qq ~ 0.5 GeV, it is shown numerically 
that these effects are large enough at Q ~ 2 GeV that they have to be taken into 
account in the next-to-leading-order analysis for exclusive processes. Alternatively, by 
introducing a new distribution amplitude that evolves more smoothly, this logarithmic 
modification can be included in the hard-scattering part of the considered process. 
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O-i. 1 Introduction 

The perturbative approach for hard exclusive quantum chromo dynamic (QCD) processes was 
developed for more than one decade [I], @> @, 3 ( see Ret- f° r reviews). In this approach, the 
scattering amplitude at large momentum transfer Q 2 factorizes as a convolution of process- 
independent distribution amplitudes, with a process-dependent perturbatively computable 
hard-scattering amplitude. By using the leading-order perturbative QCD (pQCD) analysis, 
which was performed for a large number of exclusive processes including mesons and baryons, 
the qualitatively behavior for large Q 2 could be well understood || [/J. However, using the 
asymptotic distribution amplitudes, which follow directly from the solution of the evolution 
equation, results in predicted normalizations for the elastic form factors at experimental 
accessible momentum transfer that are too small; in the case of the magnetic nucleon form 
factor, this provides the opposite sign. 

From deep inelastic scattering, where the application of pQCD is generally accepted, it 
is known that the used parton distribution functions for accessible Q 2 are far from their 
asymptotic form where all higher moments m n , i.e., n > 0, vanish. It is therefore expected 
that for the exclusive processes at accessible momentum transfer, the distribution amplitudes 
are nonasymptotical. Choosing distribution amplitudes that are enhanced in the endpoint 
region (and asymmetric for nucleons) provides the observed normalization and sign for the 
elastic form factors. 



Reference || argues that choosing such enhanced amplitudes provides inconsistencies that 
affect the importance of higher twist contributions, as well as of perturbative nonleading- 
order terms, and so the pQCD approach to elastic form factors probably is not self-consistent. 
(A second point widely discussed in the literature is the nonperturbative contribution from 
the hadronic wave function || ||.) Phenomeno logical methods, such as (1) introduce a gluon 
mass, (2) freeze the running coupling constant for small virtuality fLOf , or (3) suppress the 
endpoint region by suitable distribution amplitudes or by a cutoff [TT[|, are used to improve 
the stability of the pQCD approach. Recent incorporation of Sudakov suppression has shown 
that the pQCD approach for the pion form factor is self-consistent for a momentum transfer 
of Q ~ 20A QCD 111 (see also Ref. @). 



The validity of the pQCD approach for exclusive processes can also be studied by direct 
calculations of higher twist and perturbative nonleading contributions. It appears that higher 
twist analyses have not been achieved quantitatively. The stability of the perturbation 
theory has been investigated neglecting the evolution of the distribution amplitude by next- 
to- leading-order calculations for the pion transition form factor Ref. [14], |I5] , the pion form 
factor |Tj| [T7|, and the two-photon processes 77 — > M + M~(M = 7T,K) |T8[] . Discrepancies 
in the one-loop approximation of the hard-scattering amplitude for the pion form factor were 
clarified in Ref. [|TJ], |2D[ . The next-to-leading-order correction to the pion form factor and to 
the processes 77 — > M + M~ are rather large at accessible momentum transfer. 

Including the evolution of the distribution amplitude in these analyses requires the solu- 
tion of the differential-integral evolution equation, which can be done by using the moment 
method. The corresponding two-loop approximation of the integral kernel was computed by 
different authors and the obtained results agree with each other |2"I| . It has been confirmed 
that the computed evolution kernel is consistent with the Gribov-Lipatov-Altarelli-Parisi 
kernel |EJ and with conformal symmetry breaking in massless gauge field theories ||23[| . Be- 



cause of the complicated structure of the evolution kernel, only the first few moments of 
the evolution kernel had been computed numerically [53 j. Based on this incomplete com- 



putation, it was believed that the next-to-leading-order correction to the evolution of the 
distribution amplitude and the contribution of this correction to the pion form factor are 
rather small p4|, g5j. 

Recently, using conformal constraints, the complete formal solution of the evolution equa- 
tion in next-to-leading order could be obtained without knowing the evolution kernel by a 
one-loop calculation [2"3fl. This paper studies this solution in detail and shows that the evo- 



lution of the distribution amplitude must be included in the next-to-leading order analysis. 
Section 2 reviews to leading order the evolution equation of the distribution amplitude and 
the solution in terms of the conformal spin expansion. The evolution of the distribution 
amplitude in next-to-leading order for fixed a s is studied in Section 3. This includes a de- 
tailed investigation of the large n behavior for the next-to-leading order corrections to the 
eigenfunctions (p^(x,a s ) and eigenvalues 7 n (a s ) of the evolution kernel. Numerical results 
for the evolution of the asymptotic, the Chernyak-Zhitnitsky two-hump, and another con- 
vex distribution amplitude are presented. Section 4 analyzes the solution of the evolution 
equation in next-to-leading order with running coupling, showing by numerical computation 
that the next-to-leading-order corrections are also large in this case. Section 5 discusses the 



obtained result, comparing it with a previous result [24], and presents the conclusions. 



2 The Distribution Amplitude and their Evolution 

The distribution amplitude ip(x, Q 2 ) is the probability amplitude for finding a valence quark 
[antiquark] with light cone momentum fraction x [1 — x] in the pion probed at large mo- 
mentum square Q 2 |1|. This amplitude can be defined as expectation value of renormalized 
nonlocal light cone operators [I], |22| 

<p(x } Q 2 ) = Z" 1 /^ exp[»«(nP)(2x-l)] (0\O( K ;n)\P) | ^ 2 = Q2 , (1) 

where for simplicity, the renormalization point /j, 2 is set equal to the large momentum transfer 
Q 2 (this choice is not optimal with respect to the factorization scale setting for the considered 
processes; however, it is sufficient for the following discussions). The light ray vector n is 
chosen as n = (n + — 0, n_ = 2, 0j_) so that nP = P + , \P) denotes the pion state with 
momentum P, and 

0(n]h) = :$ d (-Kn)j 5 (nj)U(-Kn,Kh)i) u (Kn): (2) 

is the light-cone operator with the flavor content of the considered pion. The path ordered 
phase factor U(—Kn, Kn) ensures the gauge invariance of this operator. The pion decay 
constant f w = 133 MeV introduced in (1) guarantees the normalization |]J 



dx^(x,Q 2 ) = fc 1 (Q\:M0) l5 (n 1 )MQ)--\P)/nP = 1. (3) 
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Analogous to a quantum mechanical ground state, it is to be expected that <p(x, Q 2 ) can be 
chosen positive. Notice that because of charge conjugation invariance, the symmetry relation 
ip(x, Q 2 ) = <p(l — x, Q 2 ) holds true. 

The evolution equation for <p(x, Q 2 ) derived in Ref. |I|, can also be obtained in a 
straightforward manner from the renormalization group equation of the nonlocal operator 
0{K-n)^ 

Q 2 ^Lp(x,Q 2 ) = J Q 1 dyV(x,y;a s (Q 2 )) V (y,Q 2 ) , (4) 

where a s = g 2 /(An) is the QCD fine structure constant. The evolution kernel V(x, y; a s ) = 
(a s /27r) V(°'(x, y) + (a s /2ii) 2 V' 1 '^, y) + • • • has been computed perturbatively in one- and 
two-loop approximation by using the dimensional regularization in the modified minimal 
subtracted MS-scheme |21| . 



The evolution equation (|j) can be solved by conformal spin expansion 

^,Q 2 ) = t' il ^r 1 c 3 j 2 (2x-i) (o\o n (^)\py; d , (5) 



N n 



(n+l)(n + 2) 



4(2n + 3) ' 

where the sum runs only over even n [to ensure the above mentioned symmetry of (p(x, Q 2 )]- 
Here, (0\O n (Q 2 )\P) re = J dxC^ 2 (2x — l)(p(x,Q 2 ) are reduced expectation values of local 



operators that in leading order do not mix under renormalization 13, 36 1. In the free field 



theory, these operators labeled by the conformal spin form an infinite irreducible represen- 
tation of the so-called collinear conformal algebra, which is a subalgebra 0(2, 1) of the full 
conformal algebra 0(4,2) [27J. The Gegenbauer polynomials Cf/ 2 of order 3/2 form an or- 
thogonal and complete basis in the space of quadrate integrable functions with the weight 
(1 — x)x. Thus, expansion (|5|) converges if <p(x, Q 2 ) vanishes at the endpoints of the interval 
[0, 1]; see, for instance, Ref. |32|. This condition is automatically satisfied ||. 

The Q-dependence of {0\O n (Q 2 )\P) re can be determined from the evolution equation 

Q2 4v (0|O,.(Q 2 )|P) = \ £' i nk {a s {Q 2 )) (o\O k (Q 2 )\P) , (6) 

where the anomalous dimension matrix ^ n k = (a s /2iT) 7^<5 n fc + (a s /27r) 2 7„ fc +• • • is diagonal 
in one-loop order. In general, Poincare-invariance of the theory assures the triangularity of 
the matrix 7 : = {7„fc}. The eigenvalues 7„ = 7„ n are identical with the flavor nonsinglet 
anomalous dimensions known from deep inelastic scattering (moments of the Gribov-Lipatov- 
Altarelli-Parisi kernel). In leading order the solution of (|6|) is given by 

7i°V/3o 

o|o„(0-)|P) = ['—-^) (a\o n (Ql)\P) , (7) 
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A)1b(QVA 2 ) ' 
where Qq is an appropriate reference momentum, A is the QCD scale parameter, 



7-, 



1 0) = C, 



3 + 



2 n+1 - 



- 4 £ 7 



n + l)(n + 2) <-? x % 



Cp = 4/3 and (3 = (11/3) Ca — (2/3) rif, with rif is the number of active quarks and Ca = 3. 
Since 7^ < for n > 0, (0\O n (Q 2 )\P) vc decrease [see Eq. (^)] with increasing Q 2 , so 
that all harmonics with n > will also be suppressed. Furthermore, current conservation 
implies % — so that from Eq. (5) the asymptotic distribution amplitude follows: 



<p as (x) = lim <^(x,Q 2 ) = 6(1 -x)x, (9) 

Q 2 ^oo 

which does not evolve in leading order. 



In next-to-leading order the operators mix under renormalization with each other. Thus 
the evolution of (0|O n (Q 2 )|P) rc is determined by an infinite coupled first-order differential 
equation system. Since the anomalous dimension matrix is triangular, this system can be 
perturbatively solved, resulting in a behavior qualitatively different than the solution from 
leading order. 

For instance, if the initial condition is set as (p(x, Qq) = (p as (x) = 6(1 — x)x at the 
reference momentum square Q$, then all higher harmonics will also be excited. In the limit 
Q 2 — > oo, these excitations disappear, returning to ip as {x). This effect is investigated more 
generally and quantitatively in the following two sections. 

3 Next to Leading Analysis for Fixed Coupling Con- 
stant 

To see the essential features of the next-to-leading-order correction, consider first the solution 
of the evolution equation for fixed coupling constant a s . In this case, the mentioned excita- 
tion of higher harmonics by evolution will not disappear in the asymptotic limit Q 2 — > oo. 
Expansion of <p(x, Q 2 ) with respect to the eigenfunctions ip e J{x,a s ) of the evolution kernel 
V(x, y, a s ) provides immediately the solution of the evolution equation: 



<r 



X / < o2\7n("s)/2 



(x,Q 2 ) = £'^fo«.) gj (0\O n (Q 2 )\P) IC . (10) 



n=0 \^ ! 



The next-to-leading-order corrections to the evolution enters as a two-loop contribution of 
the eigenvalues j n (a s )/2 and as a s corrections to the eigenfunctions (p%f(x, a s ). The two-loop 
corrections of the eigenvalues are well known from the next-to-leading-order analysis of deep 
inelastic scattering |29|. A closed expression for the a s corrections to the eigenfunctions can 
be derived from conformal constraints and a one-loop calculation of the special conformal 
anomaly in Ref. |23| (Here, the result is re-expressed by a linear combination of Lerch 
transcendent <p[x, l,i], and taking into account the term proportional to /3 .), 

tf{x, a s ) = (-IT 2 ~§^$ £; * 1+n (1 " -) 1+n (l + g Fn(x) + 0(a 2 s )) , (11) 



where 



FJx) 



(7i 0) - A,) 



- In (x(l - x)) - ip(2 + n)+ ip(A + In) 



C\ 



In 2 



2 (1 — X 



1+n 



- E 



1 1 + 5, 
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i=l 



2 + n 



l-z, l,i) + 0(x, l,f 



'( 3 + 2n)( 7 j B + ^(2 + n)) vr 2 ' 

2 J : = : ^ + ip'(2 + n) - — 



(1 + n) (2 + 



n 



(12) 



where ip(z) = d\n(T[z])/dz, j E = 0.5772, . . . , and <j>(x, l,i) = T,%L x k /(i + k). The term 
proportional to 7^°- 1 in Eq. ( |T^ ) can be obtained directly by assuming a nontrivial fixpoint 
a*, i.e., /3(a*) = 0, from a conformal operator product expansion |JI|. I thus refer to it 
as conformal symmetry predicted part. Conformal symmetry breaking by the /3-function 
provides a shift of the anomalous dimensions j^ — > 7^°- ) — /?o- The remaining term in 
Eq. flTJD is proportional to the color factor CV, and can be interpreted as an 'additional' 
conformal symmetry breaking term that comes from the renormalization of the conformal 
operators in gauge field theory. 



3.1 Corrections to the eigenfunctions 

Consider the asymptotic limit Q 2 — > oo. As in leading order, the asymptotic distribution 
amplitude is completely determined by the eigenfunction </?q 

(p% s (x } a s ) = ip e f (x,a s ) , 



a. 



6{l-x)x\l + -^<C F 



In 2 



1 — x 
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In ul — x)xj H — 
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The term in cp^ proportional to /3o, gives a logarithmic modification. It is very interesting 
that the conformal symmetry breaking term provides an unexpected In 2 modification of the 
endpoint behavior. The a s correction to the asymptotic distribution amplitude ( |T3"D is shown 
in Figs. 2(a,b). 

We next study quantitatively the a s contributions for the eigenfunctions with arbitrary n. 
For this purpose, it is technically more convenient to deal with the following representa- 
tion |23[: 



if e J(x,a s ) 



^^ CT{2x - 1) + | tfW(x) + 0{a]) , 



x)x 



£ i^ F ^c!(2x-l)c« 



k=n+2 



N h 



kn ' 



(14) 



where 



'kn 



(2n + 3) ( 7 (°) - A) + 4A fcn ) 2(2n + 3) (A fcn - ^{k + 2) + V(l) 



(fc - n) (k + n + 3) 



(n + 1) (n + 2) 



and 



.4 



/,:?( 



cv 



V f fc + n + 4 ^| - V (^—^] +24j(k-n)-^(k + 2)- ^(1) 



(15) 



are only nonzero if fc — n even. To comprehend these a s contributions quantitatively, consider 
the amplitude at x — 0.5. However, since the In 2 ( (1 — x)/x\ term in Eq. (13) disappears at 
x = 0.5, it is clear that the large contributions of the endpoint region will be dropped out. 
Nevertheless, from Eq. (0) and C 2 3 f (0) = (-l) (n) T(3/2 + n) / (r(l + n) T(3/2)) 0, the 

relative contributions r$ = ip^ e ^(0.5) / </4°) e ^(0.5) increase logarithmically with n, and are 
of order 2 for n = 10 (/3 = 0), respectively, for n = 2 (/3 = 9). 

To take into account the missed logarithmic modification in the endpoint behavior, it is 
more reasonable to use the following quantitative measure for the 0(a s ) contribution: 



RJa s 



o (1 — x) X 



J/2 



2^ 



R { n ] + 0(al 



(16) 



Figure l(a,b) shows that this analysis provides qualitatively the same n-dependence as for 
r£\ and that the a s contributions are now larger. Moreover, the following are common 
features of the a s corrections to the eigenf unctions: 

• For n = and /?o = 0, only the 'additional' conformal symmetry breaking part gives a 
contribution, of order a s /2n. For /5 7^ 0, this term is partly cancelled. 

• Contributions from the symmetry predicted and breaking parts have different phases, 
so that the net-contribution is smaller. 



In the case of /3q = 0, the minimum is at n=6. For (3$ 7^ 0, this effect is washed out. 
For small n and (3q = 0, the corrections are small. 

The relative corrections are growing logarithmically, 

r« ~ 0.347A) - (2.71 + 1.39 ln(2 + n)) C F , 



^2 

'F 



RW ~ [0.411/3 2 + (54.7-35.9 ln(2 + n)+ 6.58 ln 2 (2 + n))C| 

+ ( - 8.98 + 3.29 ln(2 + n))/3 C F ] (V2) , 
and in the limit n — > 00, the relative corrections are independent of Pq. 



Later, the evolution of <p(x, Q 2 ) will be computed numerically. For this purpose, it is 
necessary to know how well the partial sums 



vT\*) 



1$ {l-x)x n z /2 

k=n+2 



N k 



Ck (2x — 1) c kn 



(17) 



approximate the functional series (14). This is also important for the case of running cou- 
pling, where the partial waves beyond the leading order are given by the functional series 
that have convergence properties similar to the series for the eigenfunctions. The relative 
deviation from tpf/^ (x) can be measured by 



Jo 1 dx^ 1 \xY/(x(l-x)) 



1 , ^{ (1) (*) 2 

ax 



n+2i (r^) 



"" \\ /„' dx ^«')(x) 2 / (x(l - I)) ' J» *(1 - x) t _„ +2 

Numerical computation shows that for n — 0, where /3q = 0, the deviation is 43% for 
i — 1, about 10% for i = 5, and about 1% for z = 21. In general, to get the same deviation 
for n > 0, a larger number of terms is taken into account; e.g., for n = 4 the deviation is 
50% for i = 5, 10% for i = 19, and 1% for i = 82. In Figs. l(c,d), the n-dependence of the 
deviation is shown for the cases that keep (c) two terms and (d) ten terms of the expansion 
(|17|). To remain under the 3% level for n < 500 it is necessary to keep 50 terms. The 
asymptotic expansion of A ni for large % and n, where % <n, 



A • ~ 

L - x ni — 



\ 



0.5- 



2C F [2.96 + ln(l + i)] 



/5 -C F [0.692-41n(2 + 



n. 



48.7C F - 8.42C F /3 - C F [0.692 - 41n(2 + n)] 



+ 



(19) 



0.411 + 



C F [0.692-41n(2 + n)] 



is proportional to 1/yl + i for fixed n. Furthermore, A ni increases with n and has the limit 
linin^oo A ni ~ 0.78/\/l + i. In this limit there are much larger values when n is moderately 
large; e.g., n ~ 100: 



lim A 



nl 



0.49 



lim A 



ralO 



0.235 



lim A 



n50 



0.11 



lim A 



n5000 



0.011 . (20) 



To approximate the logarithmic endpoint behavior of (p% (x, a s ), a much larger number of 
terms than suggested from the previous analysis should be taken into account. For situations 
where the endpoint behavior is crucial, e.g., for the next-to- leading-order analysis of the 
elastic pion form factor, it is better to use the following integral representation p3f : 



fn(x,a s ) = 
where 

c {l) (x,y) 



J Q 1 dy(5(x-y) + ^-c^\x,y) + 



[i-y)y 



C"J\2y-l), (21) 



(I-V)(^S(x,y) 



dzS(x,z)V^(z,y) + \g(x,y)l 



[9(x,y)]+ 



g{x,y)-8{x-y) / dz g(z,y) 
Jo 



In 1-i 
C F 9(y-x) -A >LL 



(x-y) 



x 

y 



1 — x 
1-2/ 



(22) 



Furthermore, the convolution with the kernel S(x, y) generates a shift of the Gegenbauer 
polynomial order 

d f(l -x)x) 



f 1 dyS(x, y) Q-M <#*{2y - 1) = ±. ^ ?fl C^(2x - 1) , (23) 

Jo JM n dp JS n \p=° 

and the operator V projects on the diagonal part of the expansion of a function f(x, y) with 
respect to Cf /2 ; i.e., Vf(x, y) = ££ (1 ~ x)x/N l Cf /2 (2x - l)f tl Cf /2 (2y - 1), where fa with 
< i,j < oo are the expansion coefficients. Although the operator V and the kernel S(x,y) 
are only defined implicitly, Eq. (|22|) is nevertheless helpful to convolute c^ l '(x, y) with a given 
hard scattering amplitude. 

Finally, from Eq. fl2~l|) , the a s corrections to the eigenfunctions can be written as convo- 
lution 

S e ^(x,Q 2 ) = £-£dy (c«(x,y)+ •••) V d (y,Q 2 ) , (24) 

where the partial waves of (p d (x, Q 2 ) are given as Gegenbauer polynomials 



V d M 2 ) = E' i V^ Cl'\2x-l) [%) (0\O n (Ql)\P) ■ (25) 

n=0 iV ™ \V0/ 

A further advantage of the representation (^) is that the above mentioned excitation of 
higher harmonics is now completely included in the kernel c^(x, y). 

3.2 Corrections to the eigenvalues 



The two-loop corrections to the anomalous dimensions 7„(a s ) are given in Ref. [29]. As in 
one-loop order j^ < for all n > holds true. Thus, if these two-loop corrections are 

resumed in 

{Q/Qo) (a s /2n)^ + (a s /2,y^^ 

the (modified) partial waves for n > will be more strongly suppressed than in leading 
order. The relative two-loop corrections to 7 n (a s ) are about 4.5a s /(27r) [4a s /(27r)] for all 
n > and nj = 3 [nf = 4], giving a correction of 20% for reliable values of a s ~ 0.35. The 
relative correction to the evolution of the distribution amplitude is probably of the same 
order. (This kind of correction does not appear directly in the evolution of the asymptotic 
distribution amplitude, so that in this case they are much smaller.) 

If the corrections arising from the eigenvalues are expanded with respect to a s , it is 
possible to write these corrections as convolution with the leading order solution of the 
evolution equation, 

5 ev ^(x,Q 2 ) = g) 2 ln (¥) £dy W\x,y) <f LO (y,Q 2 ) (26) 

where V d ^(x,y) = W^(x,y) is the diagonal part of V^(x, y). 



Although the kernel V^'(x,y) is known in a closed form, it seems a more difficult task 
to extract the diagonal part V d ^(x,y). A reasonable approximation can be found from the 
fact that 7W grows like 7^, i.e., only logarithmically, for increasing n. The simple form of 
the asymptotic expansion 



7: 



os(0) 



yw(i) 

In 



-5.3333 ln(2 + n)+ 0.9215 , 

-(33.237 - 2.963n f ) ln(2 + n) + 15.315 - 1.4363n/ , (27) 

and the eigenvalue equation [which is known from the one-loop approximation of V(x,y)} 

dy[v b (x,y)] + (l-y)yC 3 J 2 (2y-l) = 2(1 - lE - ^{n + 2)) (1 - x) x C 3 J 2 (2x - 1) , 



[vb(x,y)] + 

v b (x,y) = 9{y-x) 



v b (x, y) - S(x - y) / dz ^(z, y) 
Jo 



m 



X 



+ 



1 — X 



y{y-x) [y^l-y 
where ip{n + 2) = ln(n + 2) + 0(l/n) for large n, allows us to reexpress Eq. (p6|) as 

f 

/o 



cPV(x,Q 2 



it) ln Ui 



dj/ (o<J(a; -y) + b[v b {x, y)] + ) V LO (y, Q 2 ) + R(x, Q' 



,(29) 



where a ~ 0.6315 — 0.0918rij and 6 ~ 8.309 — 0.7408nf. The terms in the sum representation 
of the remainder 



R(x,Q 2 



dy (V« (x, y) - a8(x - y) - b[v b (x, y)] + ) ^(y, Q 2 



1 g, ii_^ C 3/ 2(2a; _ 1} ( 7 (1) _ T a S (l) 



n=0 



'Q_ 



2X a s7 i 0) /(47r) 



o|o„(gg)|p 



(30) 



red 



are additionally suppressed by 0(l/n). Thus, for the same accuracy, the approximation of 
R(x, Q 2 ) by a partial sum requires less terms than the approximation of (p LO (y, Q 2 ) itself. 



3.3 Complete next-to-leading order corrections 

In the asymptotic limit, each given distribution amplitude (p(x, Qq) at reference momentum 
square Ql extends into the asymptotic distribution amplitude (|13"D . Thus, in this limit, the 



relative next-to-leading-order correction 



<p NLO (x) - ip LO (x) 



f(p (x) is uniquely given by 



^asNLO^ _ (f asLO^ 

ip asLO (x) 



£^ 



In 2 



1 — x 



X 
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A) 



In ( (1 — x)x 



,(31) 



so it is large and enhanced in the endpoint region. The next-to-leading-order contribution 
by the evolution of the distribution amplitude is also important away from this asymptotic 
limit. 



It is possible to get information about the distribution amplitude at low momentum 
transfer; e.g., Q ~ 0.5 GeV, from nonperturbative methods such as sum rules |32f and 



lattice calculation |53]. However, the obtained results are inconclusive, so it is not possible 
to distinguish between the following parameterizations: 

</ s ( x ) = 6z(l-a;) , 

tp cz {x) = 30x(l - x) (l - 4x(l - xj) , (32) 

1/2 



^ co (x) = - (z(l - 

71 V 



The function y9 co (x) which was used for the next-to- leading-order analyses of the pion form 
factor in Ref. |l7j is only one example of further convex amplitudes. (For a numerical 
calculation, (p co (x) is more suitable than broader amplitudes, which had previously been 
assumed to be more realistic.) Furthermore, it is assumed that the evolution of <p(x, Q 2 ) for 
Q > 0.5 GeV can be obtained from the perturbative solution of the evolution equation. 

The evolution of <p(x, Q 2 ) is controlled by Eq. (|I0~D , where the reduced expectation values 
(0|O n (<5o)|-P) rC are computed from the nonperturbative input ip(x,Q 2 ) ), which is assumed 
to be one of the functions in Eq. fl32|) . It follows from Eqs. (|H]) and (0) up to corrections 
of order 0(a 2 ^ 



S 1 



red , „n> CX q 



0\O n (Q 2 )\Py e = m n {Ql) - -1 Y! c » m ^QD . 

m n {Ql) = f 1 dx C 3 J 2 (2x - 1) (p(x, Ql) , (33) 

Jo 

where the coefficients c ni are defined in Eq. fllBl). 

Taking into account a sufficient number of terms in the series flTH|), the distribution 
amplitude at the factorization scale for exclusive processes assumed to be Q ~ 2 GeV can 
be obtained numerically. The number of active flavors is three, and the value for the fixed 
coupling constant is a s = 0.5. The distribution amplitude was approximated by the first 100 
nontrivial terms (i — 0, 2, . . . ,200). The corresponding eigenf unctions (pf/(x,a s ) take into 
account the (102 — z/2) terms of the expansions with respect to Gegenbauer polynomials. 
The distribution amplitude <p(xj, Q 2 ) at Q = 2 GeV was then computed for different points 
Xj, j = 0,1, ... , 70 and interpolated to a smooth function. 

It can be seen in Figs. 2(b,c,d) that the relative next-to- leading-order corrections have 
the following features: 

• Independent of the shape of <p(x, Qo), the relative next-to-leading-order corrections are 
characterized by logarithmic enhancement at the endpoints caused by both corrections 
to the eigenfunctions and to the eigenvalues. 

• For partial waves with n > 0, the corrections coming from the eigenvalues are larger 
than from the eigenfunctions. However, these corrections disappear in the asymptotic 
limit. 
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• Amplitudes enhanced at the endpoints, also have larger relative next-to-leading-order 
corrections that are negative. 

Although it was possible for the chosen distribution amplitudes to compute the evolution 
in next-to-leading-order numerically, this will be a difficult task for amplitudes that are 
broader. In addition, the complete next-to-leading-order analyses for an exclusive process can 
be done more conveniently if the next-to-leading-order correction is written as a convolution, 
with the distribution amplitude (p d (x,Q 2 ) defined in Eq. (p5|), which also evolves smoothly 
in next-to-leading-order (no excitation of higher harmonics). Since ip(x,Q 2 ) = ip d (x,Q 2 ) + 
5 ef <p(x,Q 2 ) fromEq. (g§, 

(34) 



p(x,Q 2 ) = J o 1 dy(5(x-y) + ^-c( 1 \x,y) + ---^ d (y,Q 2 )- 



Again, the excitation of the higher partial waves is completely included in the convolution 
with c^(x, y). Notice that (f d (y, Ql) may be used instead of tp(y, Ql) as an initial condition. 
In fact, this corresponds to the choice of another factorization scheme for the considered 
exclusive process (redefinition of the soft and hard parts). 

The complete a s correction to the evolution of the distribution amplitude in next-to- 
leading-order can easily be obtained from ( p4|) and (|26|): 



p(x,Q 2 



V LO (x, Q 2 ) + 6 ef <f(x, Q 2 ) + 6 ev <p(x, Q 2 ) , 



(35) 



O'.s 



dy[8(x-y) + 

o V 2,ti 



c^(x,y) + ^\n(^)v d ^(x,y) 



+ 



v LU (y,Q 2 ) 



4 Next-to-Leading-Order Analysis for Running Cou- 
pling Constant 

This section discusses the solution of the evolution equation in next-to-leading-order for 



running coupling, which was derived in |23|, |3T 

OO |-1 „Q2 7, , 

<p(x } Q 2 ) = £ <p n (x } a s (Q 2 )) exp - / - ln (g(t)) (o\O n (Q 2 )\P 



n=0 



red 



(36) 



The partial waves <p n (x, a s (Q 2 )) are now Q 2 dependent nonpolynomial functions, known as 
functional series 



ip n (x, a s (Q : 
^\x,Q 2 ) 



(l-x)x * 



N n 



C£(2x-1) + 



a s (Q 2 

2vr 



,(i) 



^ } {x^) + 



k=n+2 



' (l-x)x * 



N k 



Ci(2x-l)s kn (a s (Q 2 )) c 



(i) 

kn i 



(37) 



,d) 



where c^. n are the expansion coefficients of the eigenfunction defined in Eq. (|H^) and 



$ki 



,,(a s (Q 2 



7f ~ 7f 
7f-7i 0) + /5 



1 _/^Ql)V + (7i 0) -7i 0) )//?o 
U(Q 2 W 



(38) 
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Since these partial waves satisfy the convenient initial condition 

fn(x, a s (Ql)) = ^^ C 3 k /2 (2x - 1) , (39) 

the expectation values {Q\O n {Q^)\P) can be now simpler computed as for fixed a s 

i \ red P 

(0\O n (Ql)\P) = J dx Ci(2x- 1) <p(x,Ql) . (40) 

Because of the asymptotic behavior of 7^ = — 41n(fc + 2) + - ■ -, Skn\ot. s {Q 2 )\ approaches 1 

for k 3> n and a s (Ql) > a s (Q 2 ) [see Fig. 3(a)]. Consequently, the behavior of(p$p (x, a s (Q 2 )) 

in the endpoint region is determined by cj^J; i.e., it has the same logarithmic modification 
as cpf/^(x). 

To avoid this excitation of higher harmonics (Gegenbauer polynomials) by the evolution, 
a new distribution amplitude analogous to the case for the fixed coupling constant is intro- 
duced that satisfies a diagonal evolution equation (the corresponding evolution kernel has 
to be diagonal with respect to Gegenbauer polynomials). A formal representation for this 



transformation kernel was given in [24 



W = r dtexp {-(#, - V {0) )t} <g> [(J - V)V (1) ] ® exp {-V {0) t} , (41) 

but, as was pointed out, this representation cannot be used for explicit calculations. Hope- 
fully, changing the factorization scheme for the exclusive process under consideration will 
allow us to factorize the process amplitude in terms of the desired diagonal distribution 
amplitude ip d (x,Q 2 ). 

For the numerical study of the next-to- leading-order corrections, assume that the distribu- 
tion amplitude at Qo = 0.5 GeV (A^ 3 ^ = 0.4 in next-to-leading-order; i.e., that a s (Ql) ~ 0.9) 
can be parametrized by one of the functions in Eq. fl3"2|). The amplitudes are evolved to a 
scale Q = 2 GeV, where a s (Q 2 ) ~ 0.3. The number of active flavors is three, taking into 
account the first 100 nontrivial terms in the partial sums for both series ( |3"E| ) and ( |3"7| ) [the 
asymptotic (Chernyak-Zhitnitsky) distribution amplitude requires only 1 (2) term(s) in (|36|)1. 
The result in Figs. 3(b,c,d) shows that the relative next-to-leading-order corrections for run- 
ning coupling have qualitative and quantitative features similar to those in the case of fixed 
coupling discussed in Section 3.3. 

5 Summary and Conclusion 

This paper has shown that the (relative) next-to-leading-order correction to the evolution of 
the pion distribution amplitude is rather large, especially in the endpoint region, and that in 
this region the negative corrections are larger for enhanced amplitudes. The a s correction to 
the partial waves comes from the off-diagonal matrix elements of 7^; it can be interpreted 
as excitation of higher harmonics (Gegenbauer polynomials) by evolution, and appears as 
ln(x(l — x)) and ln 2 (x/(l— x)) terms. The two-loop contribution to the anomalous dimension 
7„ is for n > much larger than the off-diagonal matrix elements of 7,^; i.e., about 20% 
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of the one-loop approximation. However, the exponentiation of the two-loop contribution 
provides a larger suppression of the corresponding harmonic's as in leading order (expansion 
with respect to a s provides a large (negative) excitation of the harmonics). 

The obtained large next-to-leading-order correction seems to contradict a previous anal- 
ysis |24]| , where it was found that this correction is rather small. The explanation for this 
discrepancy is that (1) only the first few expansion coefficients c n k were taken into account, 
and (2) the authors looked only to the evolution of <p(x, Q 2 ) at x — 0.5. Furthermore, the 
reference momentum chosen for use in Ref. [24] was Qq = 10A^ 3 ^ = 1 GeV in leading order; 
i.e., a s (Ql = lGeV 2 ) ~ 0.3. Such a choice provides a much smaller next-to-leading-order 
correction for Q = 1.25A® (because a s (Ql = 1.25 2 A 2 ) ~ 0.9, perturbation theory should 
be valid for the evolution of the distribution amplitude). Using a popular parameterization 
at lower reference momentum (e.g., Qq ~ 0.5 GeV) provides logarithmic correction, which 
should be included in the input amplitude at a higher reference momentum. 

The question of whether to include an a s suppressed logarithmic correction to the input 
amplitude <p(x, Ql) can be avoid by chosing a distribution amplitude that evolves smoothly, 
with no excitation of higher harmonics by evolution. The amplitude (p d (x, Q 2 ) satisfies an 
evolution equation where the corresponding evolution kernel V (x, y) is diagonal with respect 
to Gegenbauer polynomials. Consequently, in such a factorization scheme, the contribution 
responsible for the mentioned excitation of higher harmonics is now included as the a s 
correction to the hard scattering amplitude of the considered process. 

Because of the size of the discovered correction and its dependence upon the input ampli- 
tude, the evolution of the distribution amplitude has to be included in the next-to-leading- 
order analysis of exclusive hard momentum processes. For large enough Q 2 , the Sudakov 
suppression can be neglected so that, using the known expressions for the hard scattering 
amplitudes of the pion transition form factor and the electromagnetic form factor, it should 
be straightforward to re-analyze the next-to-leading-order corrections for these processes. 
Because of the large number of Feynman diagrams, the a s correction to the hard scattering 
amplitude for the 77 — > M + M~ processes for the case of equal momentum sharing was only 
computed numerically. It should nevertheless be possible to estimate the size of the correction 
coming from the evolution of the distribution amplitude. A general next-to-leading-order 
analysis for arbitrary distribution amplitudes requires an analytical calculation of the hard 
scattering amplitude (448 diagrams). 
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Captions 



Fig. 1. Values of r^ (boxes) and R^ (circles) are given when (a) (3q is set to 
zero, and (b) I3q = 9. The difference of the relative value of the conformal symmetry 
predicted part (upper half plane) and of the relative value of the 'additional' conformal 
symmetry breaking term (lower half plane) is shown by filled boxes (circles). The lines 
represent the corresponding asymptotic expressions. (Subasymptotic terms were also 
taken into account for the approximation of R^ .) The relative deviation of the partial 
sums (f^k (x) from the exact a s correction y?^ 1 - 1 (x) is given (c) for k = 2 and (d) for 
A; = 10. 



Fig. 2. Evolution of the pion distribution amplitude for fixed a s = 0.5 and three 
active flavors. As nonperturbative inputs, three distribution amplitudes defined in 
Eq. (32) are chosen at the reference momentum scale Qq = 0.5 GeV. They are shown in 
(a); (f as (x, Q 2 ) in leading order (solid), (p cz (x, Q 2 ) (dashed), (p co (x,Q 2 ) (dash-dotted), 
<p as (x,Q 2 ) in next-to-leading-order (dotted). The relative next-to-leading-order cor- 
rections at Q = 2 GeV are shown for Lp as (x, Q 2 ) in (b), for ip (x, Q 2 ) in (c), and for 
(p co (x, Q 2 ) in (d) showing that the endpoint behavior of the distribution amplitudes 
changed more drastically under evolution. The next-to-leading-order corrections of the 
eigenvalues are neglected for the dashed line, expanded with respect to a s for the dash- 
dotted line, and taken into account by resummation for the solid line. The correction 
in the asymptotic limit is dotted. 



Fig. 3. The evolution of the pion distribution amplitude for running a s is essentiall 
determined by the matrix valued function Sk n (a s (Q 2 )). (a) shows that Sk n \&s{Q 2 ) 
defined in Eq. (38) as a function of k/(n + 2) is nearly n independent, and for a s (Q 2 ) = 
0.5a s (Qo) if is almost of order 0(1). The relative next-to-leading-order corrections for 
ip as (x,Q 2 ) in (b), for ip cz (x, Q 2 ) in (c), and for <p co (x, Q 2 ) in (d) are comparable to 
the fixed coupling result. Here, a s (Ql) = 0.9, a s (Q 2 ) = 0.3, and three active flavors 
were chosen. The meaning of the solid, dashed, and dash-dotted lines is the same as 
for Fig. 2(b,c,d). 
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This figure "figl-l.png" is available in "png" format from: 
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